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In this paper we present an analysis of the magnetic toroidal moment and its relation to the 
various structural modes in i?3c-distorted perovskites with magnetic cations on either the perovskite 
A or B site. We evaluate the toroidal moment in the limit of localized magnetic moments and 
show that the full magnetic symmetry can be taken into account by considering small induced 
magnetic moments on the oxygen sites. Our results give a transparent picture of the possible 
coupling between magnetization, electric polarization, and toroidal moment, thereby highlighting 
the different roles played by the various structural distortions in multiferroic BiFeOs and in the 
recently discussed isostructural material FeTiOs, which has been predicted to exhibit electric field- 
induced magnetization switching. 
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The concept of magnetic toroidal moments in solids 
has recently received increased attention due to its po- 
tential relevance in the context of multiferroic mate- 
rials and magneto-electric coupling ±&&&i2. A magnetic 
toroidal moment represents a vector-like electromagnetic 
multipole moment which breaks both space and time re- 
versal symmetries simultaneously. It can be represented 
by a current flowing through a solenoid bent into a torus, 
or alternatively, by a ring-like arrangement of magnetic 
dipoles£ The toroidal moment has been proposed as the 
primary order parameter for the low-temperature phase 
transition from a ferroelectric into a simultaneously ferro- 
electric and weakly ferromagnetic, i.e. multiferroic, phase 
in boracitesi 7 - In addition, the observation of toroidal do- 
mains in LiCoPC>4 has recently been reported. - This sug- 
gests that ferrotoroidicity is a fundamental form of ferroic 
order, equivalent to ferromagnetism, ferroelectricity, and 
ferroelasticity£ 

The practical relevance of ferrotoroidic order stems 
from the fact that the presence of a magnetic toroidal 
moment also leads to the appearance of an antisymmet- 
ric magneto-electric effect) 7 -^ This is particularly inter- 
esting considering the extensive current research efforts 
aimed at finding novel multiferroic materials which ex- 
hibit strong coupling between magnetization and elec- 
tric polarization. 10 ' 11 ! 12 ' 13 As suggested in Ref. H, the 
toroidal moment concept can offer useful guidance in or- 
der to identify possible new candidate systems and to 
analyze the specific nature of the magneto-electric cou- 
pling. At the moment, however, it is not fully clear how 
this has to be done in practice. It is therefore the pur- 
pose of this work to present an instructive analysis of the 
toroidal moment for an important class of multiferroics, 
and to illustrate how such an analysis can provide insights 
into possible coupling between the various order param- 
eters. Specifically, here we evaluate the toroidal moment 
for the case of the i?3c-distorted perovskite BiFeOs and 
the recently proposed isostructural system FeTiOs (see 
Ref. EH). 

BiFeC>3 is probably the most studied multiferroic to 
date, whereas R3c FeTi03 has only recently been pro- 
posed as a material that exhibits ferroelectrically-induced 



weak ferromagnetism, and thus offers the possibility 
of electric-field controlled magnetization switching. 14 ' 15 
First principles calculations show weak ferromagnetism 
for both BiFeC-3 and R3c FeTi0 3< ii^ The net mag- 
netization in these systems is due to a slight cant- 
ing of the mainly antiferromagnetically ordered Fe 
spins. This canting is induced by the Dzyaloshinskii- 
Moriya interaction . 17,18 A small magnetization has in- 
deed been observed experimentally in thin film samples 
of BiFeOs; 19 ' 20 whereas in bulk BiFeOs this effect is can- 
celed out by the presence of an additional cycloidal ro- 
tation of the antiferromagnetic order parameter^ As 
was shown by both symmetry analysis and explicit first 
principles calculations, the weak magnetization is lin- 
early coupled to the spontaneous electric polarization in 
FeTi0 3 , but not in BiFe0 3< ii^^ The analysis of the 
toroidal moment presented in the following confirms this 
fact while in addition providing a complementary per- 
spective. 

The toroidal moment t of a system of localized mag- 
netic moments A at sites rv can be written as££& 



t = - } , n x ra, 



(1) 



As described in Ref. 4, the presence of the position vector 
in Eq. ([1]) together with the periodic boundary conditions 
encountered in bulk systems lead to a multivalucdness of 
the toroidal moment, in close analogy to the case of the 
electric polarization. 22 ' 23 As a result only differences in 
the toroidal moment (induced for example by a structural 
distortion) are well defined quantities, and the multival- 
uedness has to be taken into account when evaluating 
such toroidal moment differences. A toroidal state is 
represented by a spontaneous toroidal moment t s ^ 0, 
where t s is evaluated as the change in toroidal moment 
with respect to a non-toroidal reference configuration. 
On the other hand, a non-toroidal state corresponds to a 
"centrosymmetric" ensemble of toroidal moment values, 
but does not necessarily imply that the straightforward 
evaluation of Eq. ([T]) for one unit cell leads to t = 04 
As already discussed in Ref. 0, the toroidal moment of 
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TABLE I: Coordinates of all ions i within the rhombohedral unit cell of the R3c ABO3 structure, fi = aiSi + 0202 + 0303. 
Without loss of generality we define the origin to coincide with the position of the first A cation. 
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BiFe03 evaluated in the limit of localized magnetic mo- 
ments vanishes if one takes into account magnetic mo- 
ments only on the nominally magnetic Fe sites. This is 
due to the fact that the antiferromagnetically ordered 
Fe cation sublattice in BiFe03 represents a simple rhom- 
bohedral lattice, with inversion centers located on each 
cation site, and thus t s — 0. The same holds true for R3c 
FeTi03. The symmetry-breaking required for a nonva- 
nishing toroidal moment in these systems is due to the 
structural distortions exhibited by the oxygen network 
surrounding the magnetic cations. In the following we 
will therefore assume that small induced magnetic mo- 
ments are located on the anion sites in both BiFc03 and 
FeTi03, and we will evaluate the toroidal moment cor- 
responding to these induced magnetic moments on the 
oxygen sites. Note that if the full magnetization den- 
sity would be taken into account when evaluating the 
toroidal moment, then the full magnetic symmetry of the 
system would automatically be included in the calcula- 
tion. A formalism for calculating the toroidal moment 
directly from the quantum mechanical wavefunction has 
been suggested recently' 

In this work, we are considering perovskite-derived sys- 
tems with structural R3c symmetry' i.e. the crystal 
structure found experimentally for BiFeC>3 at ambient 
conditions" i?3c FeTi0 3 (and MnTi0 3 ) can be syn- 



b) 





FIG. 1: a) Unit cell definition used in this work. To better 
visualize the orientation of the coordinate system and rhom- 
bohedral basis vectors a;, the rhombohedral unit cell is in- 
scribed into two cubes of the underlying perovskite structure. 
Only ions within one rhombohedral unit cell are shown. The 
depicted atomic positions correspond to the undistorted case, 
b) Orientation of the three glide planes £1,2,3, and of the 12 
equivalent directions for the antiferromagnetic order parame- 
ter L (arrows) in the x-y plane. 



thesized at high pressure, and remains metastable at 
ambient conditions, even though the equilibrium crys- 
tal structure in this case is the illmenite structure (space 
group i?3)" 

We use a rhombohedral setup with lattice vectors de- 
fined as: Si = (^nf-a, |a, |cV a 2 = (— ^a, |a, |c^, 

i?3 = (0, —a, |c) (see Fig. QJi). With this choice of coor- 
dinate system, the electric polarization is oriented along 
the z direction, whereas the magnetic order parameters 
will be oriented within the x-y plane. 

The positions of all ions within the unit cell are listed 
in Table HI The oxygen anions occupy Wyckoff positions 
6b of the R3c space group. It can easily be seen that u = 
v = w = Sb = corresponds to the undistorted "ideal 
perovskite" case (but in our case with R3m symmetry, 
due to the rhombohedral distortion of the lattice vectors 
for c/a ^ 3\/2). In the following it will be convenient to 
express the oxygen coordinates in a somewhat different 
form using u = So + |e, v = So — |e + ~^4>i an d w — 

80 — \e — "73^' ^ n ^ s notation, the oxygen positions 
are: 
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It can be seen that d), e, and <5o define three distinct 
distortions of the oxygen network: So represents the dis- 
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FIG. 2: (Color online) Displacement vectors corresponding to 
4> (blue/dark grey) and e (green/light grey) for the six oxygen 
anions in the R3c unit cell, viewed along the z direction. 



placement of the oxygen anions along the polar axis rel- 
ative to the A site cations, <f> represents the counter- 
rotation of the oxygen octahedra around this axis (in- 
cluding a breathing, i.e. an overall volume change of the 
octahedra), and e represents an additional deformation 
of the octahedra, which compresses one side of the octa- 
hedron while expanding the opposing side (see Fig. 

As stated above, for the completely undistorted struc- 
ture with (5o=<5s = e = </ , = the crystallographic 
symmetry of the system is Rim. In this case (and if 
we for now do not consider magnetic order) two primi- 
tive cells are included in our unit cell definition. On the 
other hand, for (j) ^ (but otherwise e — Sq = Sb = 0) 
the resulting symmetry is i?3c, with doubled primitive 
cell compared to (j> = 0, whereas either e ^ 0, So ^ or 
8b =/= (while all other distortion parameters are zero) 
leads to polar R3m symmetry (again with two primitive 
cells contained in our unit cell definition if we neglect 
magnetic order). 

Thus, only 0^0 leads to a crystallographic unit cell 
doubling compared to the undistorted case (and unre- 
lated to the magnetic order), whereas both <5o and e 
break space inversion symmetry and therefore represent 
polar distortions. In the case of <5o this is intuitively 
clear, whereas e is perhaps not immediately recognized 
as polar. However, e does indeed destroy the inversion 
symmetry of the system, and it can easily be verified that 
e also creates an electric polarization P = Z *Afj if 
the full Born effective charge tensor (see Ref7l30l) is used 
for Z* . (Here, Af\ is the change in the position of oxygen 
ion i due to 0.) 

It has been correctly pointed out in Ref. [3l| that three 
independent parameters are required to describe the com- 
plete distortion of the oxygen network within i?3c sym- 
metry. In our notation these three parameters are e, u>, 
and So- However, in Ref. [3l| the distortion e was in- 
correctly classified as non-polar, which leads to incorrect 
conclusions about possible magneto-electric coupling in 
BiFe03, as will become clear in the following. In fact, as 
pointed out in Ref. l32l. the two displacement patterns rep- 
resented by e and So (and 5b as well) correspond to the 



same irreducible representation of the original Pm3m 
space group, i.e. they have the same symmetry proper- 
ties, and can therefore be viewed as two components of 
the same polar distortion. 

We now evaluate the toroidal moment resulting from 
induced magnetic moments on the oxygen positions listed 
in Eqs. (f2aj) - f2fj) . For this we first have to discuss the 
symmetry of the magnetically ordered state. First prin- 
ciples calculations suggest that for both BiFe03 and i?3c 
FeTiOa the preferred orientation of the Fe magnetic mo- 
ments, and thus the weak magnetization, is perpendicu- 
lar to P^*£. This results in 12 energetically equivalent 
orientations for the antiferromagnetic order parameter 
L = mpci — TOF e 2 within the x-y plane, either parallel or 
perpendicular to any of the three c-type glide planes of 
the underlying i?3c structure (see Fig. QJ>). (Here, mp e i 
and TOFe2 are the magnetic moments of the two Fe cations 
within the crystallographic unit cell.) The magnetic or- 
der thus breaks the threefold symmetry around the polar 
z axis and reduces the crystallographic i?3c symmetry to 
the magnetic symmetry groups Cc or Cc', depending on 
whether the Fe magnetic moments are parallel or per- 
pendicular to the remaining c-type glide plane. (Here, C 
indicates a base-centered monoclinic Bravais lattice and 
d a glide plane combined with time reversal.) In the fol- 
lowing we will consider the two representative cases with 
L aligned cither along the x direction (Cc' symmetry) or 
along the y direction (Cc symmetry). 

If the antiferromagnetic order parameter defined by 
the Fe magnetic moments is directed along the x direc- 
tion, i.e. the magnetic symmetry is Cc', then the y-z 
plane is a c-type glide plane combined with time reversal. 
This symmetry operation poses the following restrictions 
on the magnetic moments m, (i = 1, . . . , 6) at the oxy- 
gen sites: m 6x = -m 2x , m 6y/z = m 2y/z , m 5x = -n%i x , 
m 5y/z = m ly/z , m 4x = -m 3x , m 4y/z = m 3y / z . Evaluat- 
ing Eq. ((T|) using these relations together with the oxygen 
positions (f2aj) - f2"fl) , and considering the multivaluedness 
according to Ref. |J, results in the following x component 
of the spontaneous toroidal moment: 

' = | [V3(/){mi z - m 2z ) + 3e(mi z + m 2lS ) j . 

(3) 

Here, we have imposed the additional constraint that 
TO3 Z = -mi 2 — m 2z , to ensure that J2i m iz = 0- This 
corresponds to a decomposition of the full moment con- 
figuration into compensated and uncompensated parts to 
ensure independence of the calculated toroidal moment 
from the choice of origin (see Ref. 0) • We note that both 
t y and t z are vanishing if appropriate multiples of lattice 
vectors are added to the atomic positions of the oxygen 
anions. This means that the components of the spon- 
taneous toroidal moment along these directions are zero 
(see Ref. 0)- 

The corresponding expression for Cc symmetry, i.e. for 
orientation of the antiferromagnetic vector L along the y 
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direction, is: 

4^ = | \j> ( m iz + m 2z - 2m 3z ) - V3e (m lz - m 2z )} 

(4) 

In this case the symmetry restrictions for the oxygen 
magnetic moments are: m§ x = rri2 X , m ey / z = -~m2 y / z , 
m 5x = mix, m 5y / z = -m ly / z , m ix = m 3x , m iy / z = 
—m 3y / z , and there is also a nontrivial z component of 
the toroidal moment. However, since this component of 
the toroidal moment does not contribute to the coupling 
between P and M (see below), we only consider t y . 

It can be seen that in general the toroidal moment 
in R3c BiFeC>3 and FeTiC>3 is nonzero for both possi- 
ble magnetic symmetries, and that it is related to the 
structural distortions of the oxygen network (t s = for 
e = <f> = 0). However, the full functional dependence 
of t s on (f>, e, and So can not be seen from Eqs. ^ 
and (U]), since in general the values of the oxygen mag- 
netic moments will also depend on these structural pa- 
rameters (including also Sb)- To gain further insight 
into possible magneto-electric coupling we now consider 
the case e = 0, 4> 7^ 0> i- e - the paraelectric reference 
phase with crystallogra phi c R3c symmetry. As has been 
pointed out in Refs. |33|, [bj], andQji, the presence of a lin- 
ear magneto-electric effect in the paraelectric reference 
phase will lead to a linear coupling between the sponta- 
neous order parameters M s and P s . In contrast, a linear 
magneto-electric effect in the multiferroic phase describes 
the coupling of an additional induced component of po- 
larization or magnetization to the corresponding recipro- 
cal fields (e.g. M(E) = M s + aE). It is therefore very 
important to clearly distinguish between the presence of 
a linear magneto-electric effect in the para-phase (where 
M s = P s = 0) and in the multiferroic phase (M s ^ and 
+ 0). 

We first consider the case L \\ x. For BiFeOs this re- 
sults in magnetic C2'/c' symmetry, which leads to the 
additional constraint ni2 Z — rti\ z and thus vanishing 

toroidal moment ti C x ^ c ^ =0. For FeTi03 the resulting 
symmetry is C2/c', which requires r(i2 Z = — m± z and thus 

ts\x — V3aipmi z . The difference between BiFeOs and 
FeTiOs in this case results from the different site sym- 
metries of the magnetic Fe sites within R3c symmetry. 
In R3c BiFeC>3 the Fe cation is located on a site with in- 
version symmetry, whereas in i?3c FeTi03 the inversion 
centers are located in between the magnetic cations (see 
also Ref. [H). For the case L || y the resulting symme- 
try is C2/c (BiFe0 3 ) or C2'/c (FeTi0 3 ). The additional 
constraints on the oxygen magnetic moments are ni2 Z = 
—mi z and m 3z — for C2/c symmetry and m\ z = ni2 Z 
for C2'/c symmetry, leading to toroidal moment compo- 
nents t^ y ^ c ^ = and t s C y ^ = a<j)(mi z — m 3z ), respec- 
tively. The calculated spontaneous toroidal moments for 
the various cases are summarized in Table |IIJ 

It can be seen that in the paraelectric R3c phase (i.e. 
for e = 0) only FeTiC>3, but not BiFeC>3, has a nonva- 



TABLE II: Calculated values of the spontaneous toroidal 
moment and corresponding magnetic symmetry groups for 
BiFeOa and FeTiOa in the paraelectric R3c structure (e — 0) 
for different orientation of the antiferromagnetic order param- 
eter. 
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nishing toroidal moment, and that the toroidal moment 
in paraelectric FeTiC>3 is related to the counter-rotations 
of the oxygen octahedra represented by <f>. The presence 
of this toroidal moment causes a linear magneto-electric 
effect M = aE with a a t via the free-energy invariant 
Etpm t ■ (PxM) (see e.g. Ref.|). This means that 

once the polarization in i?3c FeTiOs becomes nonzero 
(which of course reduces the crystallographic symmetry 
to i?3c) , it will induce a weak magnetization via the linear 
magneto-electric effect, consistent with the design crite- 
ria outlined in Ref. [lj. Such "ferroelectrically-induced 
ferromagnetism" via the linear ma gne to-electric effect 
has been originally suggested in Ref. |33|. 

On the other hand, paraelectric R3c BiFeC>3 is non- 
toroidal and does not exhibit a linear magneto-electric 
effect. Therefore, the weak magnetization in BiFeC>3 
is not ferroelectrically-induced and there is no linear 
coupling between M s and P s in the multiferroic phase. 
This is also consistent with first principles calculations, 
where for BiFeOs weak ferromagnetism occurs in both 
the ferroelectric i?3c and the paraelectric RZc structures, 
whereas for FeTiC>3 it occurs only in the ferroelectric i?3c 
struct ure i 14 i 16 

Note that in the multiferroic i?3c phase both FeTiC>3 
and BiFeOs exhibit a toroidal moment (according to Eqs. 
([3]) and ((4])) and thus a linear magneto-electric effect. 
This means that an external electric field will induce 
changes in both polarization and magnetization, linear 
in the external field, but only in FcTiC>3 the correspond- 
ing spontaneous order parameters M s and P s are linearly 
coupled. Such linear coupling between P s and M s is re- 
quired to achieve full electric-field control of the weak 
magnetization. As outlined in Refs. [3 and [TH a rever- 
sal of P s in FeTiC>3 induced by an external electric field 
will result in a corresponding reversal of M s provided 
the antiferromagnetic order parameter (or equivalently 
the toroidal moment) is fixed by a large enough mag- 
netic anisotropy. On the other hand, such electric field 
controlled switching of the weak magnetization is not ex- 
pected to occur in BiFeOs. 

The evaluation of the toroidal moment presented above 
allows to clearly identify which structural modes, in com- 
bination with the antiferromagnetic order, lead to the 
appearance of a toroidal moment and a linear magneto- 
electric effect. In contrast to the antiferromagnetic order 
parameter, which generally depends only on the orienta- 
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tion of the individual magnetic moments, it follows from 
Eq. |T]) that the toroidal moment contains information 
about where the magnetic moments are located as well 
as on how they are oriented. Furthermore, the toroidal 
moment is a macroscopic multipole moment that is re- 
lated to the (magnetic) point group symmetry, whereas 
a proper symmetry analysis of antiferromagnetic order 
requires a treatment based on the full space group sym- 
metry. In particular, antiferromagnetic order is not con- 
nected to any particular macroscopic symmetry break- 
ing, i.e. all 90 magnetic point groups are compatible 
with the existence of antiferromagnetic order. On a mi- 
croscopic space group level, antiferromagnetic order of 
course always breaks time reversal symmetry. However, 
for systems where the magnetic unit cell is a multiple 
of the crystallographic unit cell, a primitive translation 
of the original nonmagnetic lattice can be combined with 
time reversal, and as a result the corresponding magnetic 
point group still contains time reversal as a symmetry 
element^ In contrast, a toroidal moment always breaks 
space and time reversal symmetries on the macroscopic 
level, i.e. the corresponding magnetic point group does 
not contain neither space inversion 1 nor time reversal V 
(whereas the combined operation V can still be a symme- 
try element) . Since all macroscopic properties of a partic- 
ular crystal are determined by its point group rather than 
space group symmetry^ the toroidal moment appears 
to be a more appropriate quantity to classify macro- 
scopic symmetry properties compared to the antiferro- 
magnetic order parameter. In particular, the toroidal 
moment is ideally suited to discuss magneto-structural 
or magneto-electric coupling. For a given structural and 
magnetic configuration the toroidal moment can be eval- 
uated straightforwardly, applying the procedure outlined 
m 

Ref.H 

Using the relation between the toroidal moment 
and the magneto-electric tensor a, this allows to correctly 
identify which quantities determine the magneto-electric 



properties of the system. The same can of course also 
be achieved by a group theoretical analysis of the sym- 
metry properties of the various structural modes and of 
the antiferromagnetic order parameter. The straightfor- 
ward evaluation of the toroidal moment should therefore 
be considered as an alternative (or complementary) way 
to discuss magneto-electric symmetry that does not nec- 
essarily require the application of group theoretical con- 
cepts. 

Finally, we point out that the toroidal moment is only 
related to the antisymmetric part of the linear magneto- 
electric tensor a, whereas the symmetric part of a is 
connected to other electromagnetic multipole moments 
(see Ref. d). However, for the present case where weak 
ferromagnetism is caused by the Dzyaloshinskii-Moriya 
interaction, the antisymmetric component related to the 
toroidal moment is indeed the crucial part of a. 

In summary, we have shown that by evaluating the 
toroidal moment in the limit of localized magnetic mo- 
ments, a clear picture of the different roles played by 
the various structural distortions for the magneto-electric 
properties in BiFeC>3 and i?3c FeTiC>3 can be achieved. 
The toroidal moment can be used to characterize the 
magneto-electric properties in antiferromagnetic systems. 
Its usefulness stems from the fact that is depends on both 
position and orientation of the magnetic moments and 
from its well-defined macroscopic symmetry properties, 
which allow to use point groups instead of space groups, 
in contrast to a discussion based on the antiferromagnetic 
order parameter. 
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